Extractable work from ensembles of quantum batteries. Entanglement helps. 
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Motivated by the recent interest in thermodynamics of micro- and mesoscopic quantum systems 
we study the maximal amount of work that can be reversibly extracted from a quantum system 
used to store temporarily energy. Guided by the notion of passivity of a quantum state we show 
that entangling unitary controls extract in general more work than independent ones. In the limit of 
large number of copies one can reach the thermodynamical bound given by the variational principle 
for free energy. 
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I. MOTIVATION 

The recent interest in models of quantum engines and 
refrigerators stimulates theoretical efforts to precisely 
formulate fundamental thermodynamical principles and 
bounds valid on the micro- and nano-scale. In principle 
these can differ from the standard ones and converge to 
them only in the limit of macroscopic systems. A sample 
of references, including both general considerations and 
particular models, is given in 

This paper is about the amount of work that can be 
extracted from a small quantum mechanical system that 
is used to store temporarily energy to transfer it from 
a production to a consumption center. To do so we are 
not coupling such a quantum battery to external thermal 
baths in order to drive thermodynamical engines but we 
address it by controlling its dynamics by external time- 
dependent fields. The battery comes with its initial state 
p and own internal Hamiltonian H . The idealized process 
of reversible energy extraction is then governed by the 
system dynamics plus some fields that are only turned 
on during a certain interval [0,t] of time. This leads to 
a time-dependent unitary dynamics of the battery. We 
now wonder about the maximal amount of work that can 
be extracted by such a process. 

It has been known for a long time that some states 
can't deliver work in this way. Such states are called 
passive , Q . The maximal amount of work extractable 
from a battery is then the surplus energy of the initial 
state with respect to the passive state a p with the same 
eigenvalues as p. 

As we are dealing with small quantum systems we may 
wonder whether using processes that entangle two iden- 
tical copies of a given battery can yield a higher energy 
extraction. More generally, what happens to a large num- 
ber of copies? 



We numerically demonstrate that the efficiency of en- 
ergy extraction grows with the number of copies. Next we 
show rigorously that the maximal amount of extractable 
energy per battery asymptotically equals the energy dif- 
ference between the initial state p of the battery and the 
energy of the Gibbs state ujj with the same entropy as p. 
We indicate how to construct in principle a unitary that 
achieves this optimal bound. 



II. GENERAL CONTEXT 

The Hilbert space H of wave functions of the battery 
is for simplicity chosen to be d-dimensional and we pick 
as standard basis for % the eigenvectors of the system 
Hamiltonian 



H = ^2 ej \ with cj+i > €j 

3=1 



(1) 
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We assume here that the energy levels are non-degenerate 
which holds for a generic Hamiltonian. 

The time-dependent fields that will be used to extract 
energy from the battery are described by V{t) — V^(t) 
where V(t) is possibly only different from zero for < 

t < T. 

The initial state of the battery is described by a density 
matrix p and the time evolution of p is obtained from the 
Liouville-von Neumann equation 

j t p(t) = -i[H + V(t),p(t)}, P (0) = p. (2) 

The work extracted by this procedure is then 

W = Tr(pH) - Tx(p(r)H) (3) 

where the state at time r is related to the initial state p 
by a unitary transformation 



p(t) = U(t)pUHt) 



(4) 
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with U(t) the time-ordered exponential of the total 
Hamiltonian H + V(t): 



U(t) = Tcxp (-i ds (H + V(aj) 



(5) 



Note that by a proper choice of controlling term V any 
unitary U can be obtained for U(t). Therefore the max- 
imal amount of extractable work (called ergotropy in 
can be defined as 



:= Tr(pH) - min Tr(Z7 p H) 



(6) 



where the minimum is taken over all unitary transforma- 
tions of H. 

Following Pusz and Woronowicz and Lenard , we 
call a state a passive if no work can be extracted from cr, 
i.e. if for all unitaries U 

Tr(crif) < Tr(UaU^H). 

The following theorem then holds: 

Theorem 1 ([2, 3]). a is passive if and only if 

d 

a = s j i-?') oi with s j+i ^ s j ( 7 ) 

In other words, cr is passive if and only if it commutes 
with the system Hamiltonian and its eigenvalues are non- 
increasing with the energy. Given p there is a unique 
passive state a p minimizing Tr U p U'H. This state is 
obtained by a unitary rotation of p denoted by U p and 
has the form 



a p = U pP U; = J2rj\j)(3\ 



(8) 



where {rj} are the eigenvalues of p arranged in non- 
increasing order: r J+ i < rj. The corresponding min- 
imal energy is X)^=i r j e j an( i the maximal amount of 
extractable work is given by 

W max := Tr(pH) - Tr(a p H) . (9) 



III. A GENERAL BOUND ON AVAILABLE 
WORK 

We obtain here a bound on lf max by comparing the 
energies of the passive state o~ p and of the canonical Gibbs 
state uj-q with the same entropy as p. Recall that the 
canonical Gibbs state at inverse temperature (3 is given 
by 



LO/3 



exp(-/3F) 



(10) 



and that its von Neumann entropy is strictly monotoni- 
cally decreasing in j3 with range [0, logd]. The von Neu- 
mann entropy S(p) of a density matrix p is 



For any given density matrix ponTL there exists therefore 
a unique inverse temperature j3 such that S(p) = S(ljj). 

The relation between p and j3 is, of course, highly non- 
linear. 

We now use the variational principle of statistical me- 
chanics that asserts that the Gibbs canonical density ma- 
trix is that which minimizes the free energy: 

Tr(pH) - T X S(p) > Tr(ufpH) - T' S(ufp). (12) 



With our choice of j3 we obtain that 

Tr(pH) > Tr(a p H) > Tr(uj-^H) 



(13) 



and hence the thermodynamical bound on the available 
work 



W max < Tr(p H) - TrfoH) . (14) 

Generally, uj-^ is different from a p as cj-g and a p or p 
have different eigenvalues. Note, however, that the two- 
dimensional case is exceptional because there is a one-to- 
one correspondence between the entropy of a qubit state 
and its ordered eigenvalues. Generally it is not true that 
a product of two independent copies of a passive state 
still is passive. There is therefore a possibility of improv- 
ing over (j9|) on the amount of extractable work per copy 
for several copies of a system. In other words, by using 
entangling unitaries, one can in principle beat ©. 
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s (p) = -Trplogp. 



(11) 



FIG. 1. Energy per copy of passive state er®n p associated to 

8> n p 

In Fig.l the energies e^ 1 per copy of the passive state 
CT®n p obtained from a product state ® n p are plotted as 
dots for n — 1, 2, . . . , 40. The lower line shows the asymp- 
totic value of e™. The system is a three level battery 
with energy levels {0, 0.579, 1} and the passive state cor- 
responding to the initial density matrix has eigenvalues 
{0.538, 0.237, 0.224}. The values e^> have been obtained 
by rearranging the eigenvalues of ® n p and the n-copy 
Hamiltonian H^ n \ see ([T5]) . The maximal additional 
work that can be extracted on top of the single copy 
extractable work using entangling unitaries is the differ- 
ence between e^ 1 ) and e^°°\ We will compute this value 
in the next section. 
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IV. ENTANGLING BATTERIES 

A state a is called completely passive if ® n a is passive 
for all n = 1, 2, . . . with respect to the sum Hamiltonian 

71 

H^=J2^ (15) 

3=1 

where Hj is the j-th independent copy of H . Thermody- 
namic equilibrium is equivalent to complete passivity: 

Theorem 2 (HQ), a is completely passive if and only 
if it is a Gibbs state. 

We now consider n independent copies of our battery 
and apply the general bound (j!4[) to estimate the maxi- 
mal amount of available work per battery 

< ax := \ {Tr((®» H<rt) - Tr(<r 8 „ p H«) } 

<Tv(pH) -Ti(ujjH). (16) 

It is our aim to show that this bound is actually asymp- 
totically achievable: 

Theorem 3. 

l™ = Tr(pH) - Tv(u^H) . (17) 

Proof. The proof is based on the idea of typical config- 
urations. Assume that p is diagonal in the eigenbasis 
of H , this can always be achieved by a suitable unitary 
rotation and let {Yj} be the eigenvalues of p arranged 
in non-decreasing order. A configuration of length n is 
an n-tuple |i) = . . . , i n ) of indices in {1,2,..., d} 

with corresponding eigenvalue r^r^ ■ ■ • r< n of ® n p. Typ- 
ical configurations will be of the type |i) where the num- 
ber of times the index k S {1,2, ... ,d} occurs lies be- 
tween (rfe — e)n and (r^ + e)n for every k. The subspace 
spanned by these typical |i) has approximately dimen- 
sion exp(nS(p)) and each such configuration corresponds 
to the average energy 

iffM|i}=(5> ej )|i>+o(e). (18) 

3=1 

Now we repeat the same construction for the product 
of n copies of the Gibbs state ujj. As S(utp) = S(p), the 
typical subspaces of ® n p and ® n 0J-p have approximately 
the same dimension. Moreover, for both product states 
the probability of finding a system outside the typical 
subspaces is o(e). We can now find a unitary on ® n T-L 
that maps one subspace into the other. This unitary is 
highly non- unique, and generally differs from the optimal 



reordering given by U^ p but nevertheless produces the 
state with the energy close to the optimal one, i.e. 

|Tr([/^(®»t/ (e)t ir (n) ) - Tr(((g) n ^) H^) \ < no(e). 

(19) 

Using ([TB")) we obtain 

Tr(®"a^#( n )) <Tr(a®n p HW) (20) 

< Tr(U {e) {® n p) U^H^), 

which combined with (Tl9l) yields the final estimation 

Ti(pH) - Tr(i^H) > > (21) 
Tr(pH) -Tr(u^H)-o(e). 

□ 

Remark that a unitary transforming ® n p into 
cannot be product and must therefore dynamically en- 
tangle the n batteries. In the numerical example of Fig. 1 
the asymptotic value e^ 00 -* exactly coincides with uj-s. 



V. CONCLUSION 

The notion of maximal reversibly extractable work for 
a quantum battery motivated by the concept of passiv- 
ity is discussed. It is applicable to full quantum models 
of micro- or mesoscopic machines where work is supplied 
or extracted by a quantum system ('quantum battery', 
'work reservoir') instead of a time-dependent perturba- 
tion of the Hamiltonian. A proper definition of work 
is important to develop a consistent thermodynamics of 
small quantum systems which is relevant in nanotechnol- 
ogy and biophysics. Generally, the extractable work is 
smaller than the thermodynamical bound computed us- 
ing variational principle for a free energy. Using entan- 
glement one can in general extract more work per battery 
from several independent copies of a battery and asymp- 
totically reach the thermodynamical bound. However, 
the optimal procedures of work extraction are generally 
difficult to implement by realistic control Hamiltonians. 
An interesting problem for future investigation is to find 
efficiency bounds when practical restrictions are imposed 
on the available control mechanisms. 
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